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ABSTRACT 


минә dere tor the design of a recursive digital notch- 
fiker WIEN one coefficient is proposed. This type of filter 
can be implemented readily with CCD (charge-coupled device). 
The effect of coefficient on notch-gain, notch-width, and 
passband-ripple is examined. The computer implementation is 
tested for second, third, and sixth order notch-filters and 


the opem coefficient 1s found based upon limited computer 


wordlength. 
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I. INTRODUCTION 


Digital filters are finding wide applications in many 
diverse areas, such as speech processing, picture process- 
ing, sonar and radar system, digital control system, etc. 
A very useful type of such filters is the notch-filter 
which attenuates highly a particular frequency component 
in the input signal while leaving nearby frequencies 
relatively unattenuated. 

A linear digital filter can be realized physically in 
terms of the three basic operations of addition, multipli- 
cation and unit delay. The advantages gained by realizing 
a digital filter rather than its continuous counterpart 
lie principally in its stability, accuracy, reliability 
and flexibility [Refs. 7 and 20]. Furthermore, its 
frequency characteristics can be varied without altering 
the hardware simply by chanping the stored set of filter 
coefficients. The disadvantages of digital approach lie 
in two extra design considerations. Firstly, the choice 
of wordlength of the computer and analog-to-digital con- 
verter. These problems do not arise in the corresponding 
continuous design. A design filter operates on an input- 
sampled signal by means of a computational algorithm. It 
can be simulated on a general purpose computer or can be 
constructed with special purpose digital hardware, such 


as CCD (charge-coupled device) implementation; this part 





is beyond the the scope of research. Due to computer 
limited wordlength, the round-off error in computer program 
is produced. This can cause the filter characteristics to 
deviate considerably from the ideal. 

This research attempts to present a survey and 
discussion of the effect of the finite wordlength on the 
accuracy of digital notch-filter. The design procedure of 
digital notch-filter was to propose rational functions of 
a discrete variable which, when transformed into the fre- 
cena domain, have desired filtering characteristics. The 
design was proceeded directly in terms of the variable, z, 
associated with discrete time intervals rather than 
through the intermediate variables associated with con- 
tinuous time. 

The digital notch-filter with one coefficient was 
investigated and tested by implementing on computer IBM 


360/70 installed at the Naval Postgraduate School. 





( ۳ HE CHARACTERISTIC -TRANSFER -FUNCTION 
OFFDIGITAL NOTCH-FILTER 


A. TRANSFER FUNCTION 

As the first step in the design of the digital notch- 
filter, one often uses the bilinear transformation to 
approximate the analog filters in discrete time. This 
transformation process is perhaps the most utilized approach 
relating the continuous and discrete domains. The approach 
in this paper differs by starting directly in the discrete 
domain without the bilinear transformation. 

The approach is to uniformly sample the signal u(t) every 
T seconds, and use the resulting sequence as the input to a 
recursive digital filter governed by the general difference 


equation 


v(K) = bo UK) + by mür... + ۳۳ u (K-m) 
- ај v(K-1) - a, v(K-2) - ... - an v (K-m) (2-1) 


20 07405) denocte the values of the filter's input 
and output signals, respectively, at the K-th instant. 

For minimum noise generating, cost, and simplicity, one 
can reduce equation (2-1) to the minimum number of coeffi- 


mlentsssuchew (2-1) is still in the recursive form. 


v(K) = u(K) + M u(K-m) - a v (K-m) (2-2) 





The above simple recursive equation of the recursive 


digital filter is derived from (2-1) by letting 


bo zm 
by = b, FU 0 
and a, = a, = ... 7 84 4 7 0 


Then its corresponding transfer function has the form 


15350 zn 
m 


ии (2-5) 


а Z 
1 ^ 


һе 7/1 cn filter One wishes the filter characterized by 
(2-53) to effectively transmit all sinusoidal inputs with the 
exception of the Ou radiant/sec. sinusoid which it should 
eonpletely reject. Unfortunately, it is impossible to 


implement an ideal notch-filter 


Gain 


o 


۲۰٠٠٠٢۷٣٠٠٦٠٠٦٦٠ Gain frequency behavior of 


an ideal notch filter. 
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B. GAIN FACTOR OF A DIGITAL FILTER 
The poles and zeros of the transfer function (2-3) are 


computed as follows 


۱ T 
hence z = MD = SR don 


HSU, lG yes 
T 
z = Ma = Ma ef nth) 


m m 


The zeros and poles are respectively given by 


> T 
A 11072 (2-4) 
d m 
and 
e T 
p; 7 "Va- 57 (2-5) 
jd MESE 


From (2-4) and (2-5) the zeros and poles are all complex 
conjugate for m even and for m odd, one of the zeros and one 
of the poles will be negative real. 

Ha emots sand pole's pattern iS plotted in 


Figure (2-2). 
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Zug 
Z i) 
2 J R 
e 
S & 
zer zi 
5501160 50) Pole zero pattern. 
where 
21 is complex conjugate of 21 و‎ 
1 is sampling frequency, 
wo :x9pch-tfnseguency which the one in which the first 


zeromand pole occur. 





For m=3 it is given by 


Py = a. e 


After multiplying the numerator and the denominator of (2-3) 
by z^ and factoring, the following zero-pole format is 


obtained 


(z-2,) (z-2, ) ә (2-2,,) (2-25, ) 


H(z) = ——  ——. - Ys=s - (2.03) 
(2-рј) (z-p) ) ... CzP54,) (z-p5, ) 
where m is even and given by 
m 
ط8‎ 9۰ N 
If m is odd, called m = 22+1, the format becomes 
— ıı uu,” ) 
222 2-2 ин (|2=2 zz 242 
四 = 1 1 22 27 2 1 50 
(z-p4) (2-р| Ae (2-р,,) (2-Р,, ) (2-Р2 41) 
иһеге 22504] and Daat are negative real. 
The condition for the system being stable is: 
co z T 
> İH 5 ) | € co (2-7) 


w=0 
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and the steady-state response of any stable system with 
Jusi er funetron (2-6) to the sinusoidal input, 


u(K) = Sin KwT, is given by 


y(K) = |H ии Sin Kal 


where 


E | 
H (e^) - H (z)1 , 2 çiəT 


and 


6 = angle {H 27 


The term İH (eJ 913] is called the system gain factor since 
ara cterlzes the sinusoidal gain response property of 
system (2-6). Hence, the steady state sinusoidal response 
is readily obtained by evaluating the system transfer func- 
tion at LN where w is the radian frequency of the input 
"usoid and T is the sampling period. One can determine 

the system's gain and phase angle factor in a straight- 
forward manner by using a geometrical approach. This approach 
gives the advantage of relationship between the locations 

Of zeros and poles of the system transfer function and its 
gaum-taemewEbeNEUNOr. In the design of digital filter, it 
will enable us to design a system whose gain-factor approxi- 
mates any desired sinusoidal discrimination behavior. 


Substituting sell into O) 
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50771 E 
NS UE To UNC 


Frequency magnitude plot for m=3. 
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az lulu 01 
ə 


Figure (2-5) 


925 


Ј22 


i 
سح‎ 
Сә 


910 


Qo» 





۱ ۱ | | ۱ E 
(е97.2 )(е9-2 )...(279Т-2 ) (ed UT. 2, 7) (eJ 9^ -22 41) 


er dw er 
(e - 0 ( 6 XE 177 - و‎ ( 6 “Boş ys "Р2 241) 


(2-8) 
where 
۱ T 
2. = m” oJ (2n+1) m 
1. lu. ҹил. 
n= Onal, 27 
mue cL 
m 
22031 a : 
Similarly 
: T 
p; = var EDS 
na 
P224] Am 
let 
۱ T 
5 - 1س زى‎ mA 一 еј (2n*1)- 
1 77 
1 T 
B. ۳ ¿JwT _ nya 5217: 
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Mc E c 


В2 931 m 


Hence (2-8) becomes 


* * 
la, lla, ۰.۰ امرها‎ da, laan, 


0+۸ | (2-9) 


I&; | I8, E... 1821 1821 1821 


It is obviously the system gain factor at radian 
frequency o is seen to depend on the distance by which the 
System transfer function's zeros and poles are located from 
me Critical point esə. Thus, if one wishes to completely 
ıe lect the sinusoid sin Kol, it is necessary that the trans- 
fer function have the radian frequency w increases from 


zero to infinitive, one sees that the zero of transfer 


function (2-3) will fall at the first of the zeros of 


[(2n*1)-] radians for n=0 and repeats for n=1,2,... as in 
Pis. (2-5). 
Gain 
0 Dos Su ماد‎ 0) 


Figure (2-3) Gain frequency behavior of 


1٦٦3٦4٦٦167 21172131677 
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3 3 9 111۵ ۳016 - 261۳0 ۳66۳1۲ in Fig. (2-2) and the 
equation (2-9) one wishes to reject completely at 


: 1 
P m (2n*l)r. 


Z the distance a, must be equal to zero. 
One knows that 
m j (2n+11) 二 
Bez za b. € m 
1 m 
hence 


A. 


: T ۱ T 
;= GE n/p- eJ C2n*Ug - 


From the above equation, 


a. = 0 implies ps -] 


and for the gain of (2-9) is constant for all frequencies 


: T 
execpt one at Һи. . One should choose the coeffi- 


elont, a, as close to unity as possible; the closer an is 
comunity the closer the notch-filter approaches the ideal 
notch rıltesseae Shown in Fig. (2-1). 

e iris. (2-2), the parameter A = 1 - و سح‎ which 
controls the notch-bandwidth, when = close to one, 
A close to zero and so does the notch-bandwidth. Ideally, 
it is desired to make A as small as possible, but in prac- 
tice considerations one is not able to make ^ arbitrarily 
small ew will be discussed in section D of this chapter). 


The above considerations can be made clear by analyzing 


the transfer function (2-3). 
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° 


۲ heey ols OF MAGNITUDE RESPONSE 
It is discussed here the magnitude response as a 
function of frequency and the coefficients a and DE 


The equation (2-3) is rewritten as follows: 


21 * b 
_ 11 
H(z) = m unu (2-10) 
z us 


joT and substitute into (2-10), it is given by 


Let z=e 
- eJ oT + Gr 
e + a 
m 
Using eJ* = cosx + jsinx, one gets 


(Cos muT + b ) + j Sin moT 
A 5. (252) 
(Cos moT + a) + Sin moT 


The squared magnitude of MES will be 


7ھ 2 

۱ (Cos mwT + b + sin "mal 

IB (۱ - 7 

(Cos moT + a + sın mwT 

or 
2 
Ha. b * 2b Cos moT 

|u(e) 91) |? - 2 (2-13) 


1 an * 2a Cos moT 


Take derivative of (2-13) with respect to w and let the 


۱۳/۱۱۱۲ ۶ ۵ equal "to zero 
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۱ 2 
2 2101 ۳ ۱ [1 + dee Ad Cos moT] 


2 


由 2ma, T sin moT [1 + DE b 20 Cos mwT] = 0 


or 


(- 2mb, T) (1*a, ^) sin moT - 2ma T (1*b, ^) sin mwT = 0 


Hence 
2 2 


sin moT [a (1 * Би ) - pür ` an )] = 0 (2-14) 
100905 لا‎ term in the bracket is not equal to zero. 
The frequency at which the squared magnitude is maximum or 


minimum will satisfy the equation (2-14) 


sln muwT = 0 implies wT = 一 一 


l 051527... 
At n=0, one has wT=0, the magnitude has the maximum 


value and at n=l, نا‎ , the magnitude has the minimum value 


and wee 1s defined ds a notch-frequency; when n=2,3,.... 
the magnitude function just repeats as in Fig. (2-4). 


If one wishes to get the constant gain except at frequency 
T ۳ 

WEN > m= 155 5 a Sallie Flt (2-5)4 the coetricients an 

and bu are selected such that this result is achieved. 


To do that in equation (2-14), one lets the factor in the 


bracket always equal zero, i.e., 


20 





Gain 


WT 


0 T 27 ən 
m m m 
Figure (2-4) Gain frequency behavior with some 
values of an and Da: 
7 5 
an Hi ) - ba (iram )=0 
or 
a a a .جم‎ 
m mm m m 
153157 - а Ы. (a,b) (2-15) 
Hence a = 1 
m 


We know that the system is stable when all the roots 
of characteristic polynomial or transfer function of 
system are inside the unit circle of z-plane. 


1520000 condition, the coefficient a must be 


a <] 
m 


AS Stated in Section B, the coefficient p. has the value of 


ur ompleseliy Tejecting at the notch-frequency. One 


21 





Sees chat with HH and a > 700 qalı? (216 cannot 
be satisfied. But ve can approximate (2-15) by letting the 
coefficient 175755 close to one. 

One sees that both statements by analyzing the transfer 
function (2-3) and considering the pole-zero pattern in 
me. (222) are appropriate. 


The transfer function (2-3) now becomes: 


=m 
21100 +7 


H(z) = 
1 237 b 


(2210) 


where the subscript m is deleted and its squared magnitude 


1S given by: 


8 ۹۶ = 2 + wos mu T (2-16a) 
би 0” COS MOL 


The plot of He ey for particular value of a vs radian 
frequency w for m=3 in Fig. (2-5). The half power frequen- 


Cies are that at which the gain 1S 3dB down or 


|1) (| = $ 


or 


2 + 2 COS mwl 


> (2-17) 
1] + 4 + 28" 065 Mol 


гәјкә 


Solve (2-17) for о 


4+ 4 Cos mwT = 1 + a 2 + 2a Cos muT 


(4-2a) Cos mwT = a” - 3 
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Hence 


Cos > (2-18) 


There are CwO solutions wT and w 


1 əl vhich are symmetric 


to wT=— , the notch-frequency of filter, as illustrated in 


Fig. (2-5) 





Figure (2-6) 


where B is a notch-bandwidth. 


The notch-bandwidth B will be 


2 
Nor gy 


4 - 2a 


] 


E 1 ۲ 
2 E m Cos 


OT 


2 
pn [T - Сов 7 20 (210) 


DES 





If one inserts the parameter A=l-a into equation (2-19), 
the bandwidth will be a function of A 


2 
4T _ =1 (1 = 4A) -3 
ee 5 1 20 - ml 


Or 


2 


m _ -1 A” - 2A - 2 T 
B = [m Cos IA (2 20) 


From (2-19) and (2-20), the limit of the bandwidth B 
when the coefficient, a, approaches one (or A approaches 
2 777 15 Z€T9o0. This is valid in theory; in fact, there is 
Some IIMIt for thé coefficient, a, (also for A) in practical 


view (this phenomenon will be seen in Chapter V). 


D. LOWER BOUND OF A DUE TO FINITE PRECISION ARITHMETIC 
۳ ۰ Of the limitations caused by finite wordlength 
in all current digital computers, a real number must be 
represented by using a finite number of bits in either the 
Sid Dont or floating-point form. The error introduced 
by such a representation is discussed in this section. Only 
binary arithmetic will be discussed here. 
In a fixed-point representation, assume a number x has 
been normalized such that |x| < 1 can be represented in twos 


complement form as 


where 


LI or 0 and bo 7.۲5651656 bit. 


24 





To approximate x by a t-bit register length, one of two 


processes is generally utilized: rounding or chopping. 

In rounding, a 1 or a 0 is first added to the t-th bit, 
(Х), 1 according to whether the (t+1)-th bit, (Х),. 15 
l or 0. Then only the first t bits of the result are kept. 
i chopping, those bits beyond the most significant t bits 
are simply dropped. Since the error introduced by chopping 
is larger than that introduced by rounding, chopping arith- 
metic is not commonly used. Therefore, the discussion 
below is limited to the rounding arithmetic process. 

Let (X). be the t-bit representation of the number x. 


The error which will be introduced is defined by 


e = (X). - X 


t 


where Idco (2-21a) 


It is known that if two t-bit fixed-point numbers are 
added, their sum still contains t bits, provided that there 
is not overflow. Therefore, under this assumption of no 
overflow, fixed-point addition causes no error. On the 
other hand, the product of two t-bit numbers may have more 
than t-bits. Thus, rounding is necessary if only t bits 
are to be kept for future processing. 


In the case of a floating-point number, x is represented 


x = £2° 


25 





where f, a fraction ۱ i and 1, is called the mantissa 
and e is a binary integer called the exponent. Let x denote 
the t-bit mantissa floating-point approximation of the 
number x. The difference between x and £1(x)-(x), 1s called 
٦۳٦" ۶۱۱0714-0 error and depends on the size of x. It 5 


therefore the best measure relative to x. Then it is clear 


Elat 


(x), = x(1*e) ۳2 


where the relative error e is bounded by -2 ^ < e < 2 * 
Unlike the fixed-point case, both addition and multiplica- 
Pon an Floatıng point can introduce round-off errer. 

From the preceding section, it is apparent that a 


digital notch filter with notch at w zi can be implemented 


0 
by means of a filter vhose transfer function has a zero 
located at “7. İn order that the gain of such a filter 
in the frequency domain is near one for frequencies suffi- 
erently close to 00， it is necessary that there exist 
properly located poles which tend to cancel the effects 
of these zeros. This can be accomplished by selecting an 
appropriate A; from equation (2-20) it is obvious that the 
ls O1 lie filter 1S directly dependent on the parameter 
A. Ideally, it is desirable that A be made as small as 
possible in order to make the notch width small. 

There are practical considerations, however, which when 


taken into account, will not enable the filter designer to 


make A arbitrarily small (see Chapter V). In the design 
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01380 neech-filter of (2-16), one wishes the gain factor 


at w= to be zero. For the critical point ul 


т sufficiently 


. TT 


outside the circle of center z, F € M and of radius N 
5217077000) che gain factor for all practical purposes is 


one. 


Center 


krsuro O) A pair of first pole and zero, 
This follows since, for 20.701 outside this circle, by 
selecting A small enough, the ratio of the product of the 
distances from zeros and poles to the critical point Ex 


7> 5566۲7۰01 a constant. Since all digital computers 


27 





are limited by finite wordlengths, the desired coefficients 

will be truncated and the coefficient inaccuracies will 

result as stated at the beginning of the D of Chapter II. 
The characteristic difference equation of a digital 


filter (2-16) is given by 
v(K) » u(K) * u(K-m) - av(K-m) (2-22) 


When one wishes to implement (2-22) on a digital computer, 
the coefficient, a, will be stored with error due to trunca- 
tion as (ate), where e 15 the round-off error in (2-21a,b). 
This will rise to a shift in the actual pole locations of 
the transfer function. 


For system stability, the condition placed on (ate) is: 


a + g 1 
Or l-a > € 


hence E (2-23a) 


If one uses a t-bit wordlength, (2-23) will become: 


E: (2-23b) 


In the implementation of digital notch-filter (2-22) 
with t-bit wordlength, one must select the parameter 
A > 2° to guarantee that the notch-width is sufficiontly 
narrow and to avoid the possibility of unbounded system 
growth. It is noted that if one counts on guard bit, the 


wordlength will be (t*1) bits. Given notch-width B, one 


can calculate A from (2-20) and therefore t. 
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IDEE PACESBEPRESENTATTONS AND REALIZATIONS 


The transfer function (2-16) of an m-th order digital 
filter is the relationship between one input and one output. 
Only state and output equations with single input and single 
output will be considered. In the time invariant discrete 
System, the transfer function (2-16) can be described by 


a state variable representation of the form 
201517 x — BonÜK) (3-la) 
UK) = C x(K) 人 人) (° lb) 


where 
AS X m matrix, 
Bis anm X ] matrix 
6577: PX m matrix 
Sis] i matrix or scalar. 


Dime the TG ۲1 119 ۲ 0۶1 of (3-la) and solving for X(z) vielids 


lc 


350151 (2) 15) 


where I is an m-dimension identity matrix. Then, taking 
the z-transform of equation (3-1b) and substituting (3-2) 


for X(z) yields 


TEQ” رب‎ B + 6) U(z) 


(2) 


H(z) U(z) (3-3) 


29 





where He) = C [I = —.. B + ó (3-4) 


H(z) is the transfer function of system (3-1). The 
discrete system is considered, as stated above, for only 
one input and one output. Hence, the matrix parameters 


777) C and ó have the forms: 


811 AS ајы 
m x m 
A = 
anı Am mn 
1 
0 
mx 1 (3-5) 
B = 
0 
2۸ 7 
C = [C, 20.5 x Ca) 
6 = d 
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There are many ways to realize the digital filter (3-1). 
Some of them are the canonic realizations. The term 
"canonic" is used in the sense of minimizing the number of 
operations required. This minimization implies a minimum 
number of coefficients of the difference equations, which 
is the source of noise. 

If one wants to describe the system (3-1) by a matrix 


form as: 


S= |---+—--. (55) 


where S is called characteristic matrix system. 
The matrix parameters A, B, C and $ are respectively 


given as follows: 


一 


I 


نتم © 


B = | (5-6) 
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C= [— 0... 0] 


Cə 
1 
ده‎ | ә 


To describe a system in the matrix form, its state 


equations will have the form: 


X1 (Kri) = a E (Kİ 0 (K) 
X» (K+1) = x] (K) 
5 
Xü (K+1) = Х 1 (К) 
v (0 = (£ x, (Q) + Lu (K) 


The corresponding realization of (3-7) has the form 


5 2 








Figure (3-1) Non-canonic realization of m-th order 


digital notch-filter. 


The realization in Fig. (3-1) is not canonic; the purpose 
of this paper is to use only one coefficient. There are many 
disadvantages in the above realization. Among these are 
the production of more noise, higher cost and greater network 
complexity. Hence, the matrix. form (5-7) is a poor way to 
describe the system. This leads to the consideration of 
a search technique for the canonic realization as stated 
in section A. From this concept, one can realize the digital 


notch-filter (2-6) as follows: 


5 








Figure (3-2) Canonic realization of digital 


notch-filter. 


and the corresponding state equations of realization (3-1) 


are listed below 


X1 (K) SE (K*1) * u(K) 
X1 (K+1) = X1 (K) 
xə (K+1) = X1 (K*1) 


(5-8) 


Zo (K+1) ә dl (K*1) 


y (K) Xn (rl) mer x) (K) 


One sees the realization in Fig. (3-2); with this 
technique there is only one coefficient and it provides 
the minimum number of the coefficients of the m-th order 
recursive digital filter one can achieve, Therefore, in 
noise analysis, if one consider only the source of noise 
produced by the multiplications of the coefficients, it 


is not necessary to assume the sources of noise are 
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uncorrelated, since there is only one source of noise in 
the realization. By this advantage, in the design of high 
` order of digital notch-filter it is not necessary to imple- 
ment by cascading first or second order sections as is 
often done at present. 

The "dual" canonic realization of Fig. (3-2) can be 
derived. The so-called "transposed realization" has been 


shown by Jackson [Ref. 4] to be given in the form 


m v(K) 


u(K) X4 (K) x4 (K*1) х (К41) 


: ^ 9 


Figure (3-3) Transposed canonic realization 


9 ۰19 ۲ 11۱۵۳ mToteh-fffTterT. 


The corresponding equations are given by: 


X1 (K) = -a u(K) + u(K) 
X1 (K+1) = X1 (K) 
xə (K+1) = x] (K+] ) 
00 
X. (K+1) = Х 1 (KEL) 
KACK) ки 1) 11) 
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In the state ons (3-8) and (3-9), the value of 
the coefficient, a, is taken between Ü and 1. It is known 
that the notch-gain and notch-width of a digital notch- 
filter are functions of the coefficient a. The question 
which arises here is: what is the optimum value of the 
coefficient a which satisfies the desired characteristics 
of the notch-filters and how large is the ripple in their 
passband? 

These questions are investigated in the following 
chapters by simulating the digital notch-filter with the 
digital computer IBM 360/70. One can see the difference 
between theory and practice in Chapter III. 

For example, the gain in equation (2-9) at notch 
frequency, in theory, is a decreasing function when the 
coefficient a 1s increasing. However, in practice, some 
values of coefficient a make the notch-gain change non- 
monotonically because of round-off noise in iteration 


computations in the program. 
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IV. IMPLEMENTATION 


A. COMPUTER IMPLEMENTATION 

Ineens chapter, a digital notch-filter which is 
characterized by the state equations (3-8) is implemented 
on the IBM 360/70. Plots of the gain-frequency response 
and the time response are provided. 

The input signal is a sine-wave, i.e., a signal of 
single frequency, which sampled every T seconds. The 
corresponding gain for each input sine-wave is computed 
when the steady state is reached, by the ratio of the value 
of the peak input sine-wave divided by the value of peak 
output response. Then the input sine-wave is scanned over 
the frequency range of the notch-filter and the gain- 
frequency response is plotted. It is noted that the notch 


of the filter falls at the first zero of filter; that is, 
. TT 


at 21 = e 7 If one desires the notch-filter at و لا‎ 


radians, the sampling time T is computed as follows 


Or 


Hence 





Or T = amt. (4-1) 


ü 
where Еј = = is the notch-frequency in hertz. 


The following sections contain implementations of 
second, third, and sixth order digital notch-filters on 
EE TEM 500 computer using single precision arithmetic. 

1. Second Order Digital Notch-Filter 

The results of m-th order digital notch-filter 
in Chapters II and III are applied. The second order 
digital notch-filter is characterized by the transfer 


function 


-2 
1 + 
H(z) = ——+—, (4-2) 
Laz 


ANG ic realızed as follovs: 


d 


x5 (K*1) v (K) 









Figure (4-1) Canonic realization of second order 


digital notch-filter. 
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The corresponding state equations are given by: 


1 (K) = -a X, (RIJ TUCK) 
1 (K+1) = Xi (K) 

(4-3) 
X, (K*1) = x] (K+1) 


۷ (K) > (K+1) vi Së (K) 


If one wishes a notch frequency at fo = 60 hz, the 

equation (4-1) gives the sampling frequency T as follows 
ğı 7 1666666857 x 10” 

The figures from (4-2) to (4-10) illustrate the 
gain frequency responses in dB vith different coefficients 
of one thousand iterations. It is noted that at one thousand 
iterations the steady state of frequency response is not 
yet reached. These plots merely depict the behavior of 
the gain-frequency response with changing coefficients. 

The result of the simulating program gives the steady 
state of frequency response at about six thousand iterations. 
17 are shown in Figs. (4-2) to (4-17) corresponding to 
some typical coefficients. And the notch-gains, notch- 
widths, and maximum passband-ripple in dB of the steady 
State gain-frequency response are tabulated in Table (4-1). 
۳ 11] shows the plot of notch-width vs. coefficient 


in theory by evaluating the equation (2-19). 
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a 


The time response Of a second order digital notch- 
filter when the input is a sinusoidal signal of frequency 
60 Hz, i.e., at the notch frequency given below. In theory, 
the notch-filter with 60 Hz notch will immediately reject 
the above input signal, but in practice the time response 
sinusoidal signal will be damped in some finite time, small 
or large according to the coefficient and order of notch- 
filters one utilizes. Figures (4-18) to (4-21) show the 
time response of second order 60 Hz notch-filters with 60 Hz 


sinusoidal input signal applied. 
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- 86.25 dB 


Frequency response of second order 


digital notch-filter. a-0.90714559 
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Figure (4-2) 
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Figure (4-3) Second order 


digital notch-filter. 
az0.969901345 
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digital notch-filter. 
a=0.981 


۳۴ ھ2 


CALE::2.00£+01 UNITS INCH. 


ир 
Y-SCBLE-1.00t«01 UNITS INCH. 


45 


-032 


EU 


LS 


-096 


-097 


3 


=98 





-025 


0065 


-003 


- 80.7 dB 


UNITED INCH.‏ ۱۱ ۰ 9 8ت 
۳+ ۹ ۱۹ ۷"۹ؤ)ؤ ‏ 8 ۰ 


1500769007 Second Order digital notch-filter. 
X 


a-0.9815 


46 








dB‏ ۶٭ 


BUE ۸, " ٣ 
QGE+0L UNITS INGH. 


۶۰۶1 40 ٣ order digital 77066 77٦ 


a=0.982 


4,7 


ur 


-С01 








和 
ENT oe ٣۳ 


ıdı 


nes 
> 


1505005 179) 7 Second order digital notch-filter, 


a=0.98298407 


48 


i 


-2.2 


-t) a 


-V75 


-319 





——P — W‏ سیت ہے ھللا 


NE Se e 
x-8-——.—.-—-—. 


. 626 dE 


Second order digital 1766-0-7٠ 


a=0.99845153 


#9 


Pıiguer (4-10) 


D - 


(C1 


, d 


“x 


827 





: - 
un. 2 a ss. 


COEFFICIENDA 


Plot of notch-width vs coefficient, 


J 


s. 

hos 3-4 +. 
— -一 | 
N U7 
sece — — 
SE o üə لے‎ 
جم جم‎ 
— e 


in theory,of equation (2-19). 


50 


ZH Ut HIGIM-HO.ION 


۳ ۳ M) 


xə 


væ 


13 








— —— OS — P  — — - - — — —À— 
J A мы ) ead 412 
Ен 
2, 


Figure (4-12) Steady-state 
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It is noted that the passband ripple is large for 
frequencies very close to the notch-frequency. Further, 
about one hertz from the notch, all passband starts less 
than 0.09 dB in magnitude and smoother for going farther 
free WO t Ci. 

2 Third Order Digital Notch-Filter 
The third order digital notch-filter is charac- 


terized by the transfer function 


تی :25 


l + az 


H(z) (4-4) 


5 


and realized as follows 





У(К) 





x; (K*1) 
u (K) 


Figure (4-22)  Canonic realization of third-order 


duzıcananotcaurırlter- 
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The corresponding state equations are derived as 


CNN ME (1) + u (K) 

хј (K+1) = x4 (X) 

X, (K*1) = x, (K*1) (4-5) 
NCC CON) 

v (K) = x, (K+1) + x, (K) 


and the sampling time required for 60 Hz-notch is given as 


ll E یت رو‎ 


Similar to the results of the second order case, 
the frequency responses of a third order digital notch- 
la are Shown in Figs. (4-23) to (4-31) with one thousand 
iterations. The steady state of third order is reached in 
about ten thousand iterations. The steady-state gain 
Bz cz neyaresponses are illustrated in Figs. (4-32) through 
(4-37). 

The time responses of third order digital notch- 
ater are shonn in Figs. (4-38) to (4-41). Similarly, 
Table (4-1) of second order, Table (4-2) of third order 


is tabulated. 
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Figure (4-23) Frequency 
response of third ms 
Noten filten: 
a=0.90714359 


= 70.4 dB 


ul: “00:01 Uh] TS ۰ 
"ھ۰‎ UNITS INCH. 


64 


-002 


y 
J 


=L 


A5 


-007 





= v 
222 022 ۹9 2) | 


> 61.7 dB 


NG 
+۹۰۰) + ۰۷۳ 


Figure (4-24) Frequency response of third order 
notch-filtre. 
a=0.96901345 
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Figure (4-28) Frequency response of third order 
notch filter. 
a=0.9815 
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Figure (4-33) Steady-state frequency response of 
third order notch-filter. a=0.985 
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Figure (4-35) Steady-state frequency response of 
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Time response of third order 
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3. Sixth Order Digital Notch-Filter 

The sixth order digital notch-filter is 

characterized by the transfer function 
H(z) = (4-6) 
0 06ھ‎ 

One often realizes the sixth order digital filter 
by cascading three sections of second order filters. This 
is done to reduce the sensitivity of pole positions due 
to truncation because of round-off error of iteration 
computations in the simulation program. 

The sixth order digital notch-filter characterized 
by transfer function (4-6) is discussed here. There is 
only one source of error; the noise produced in the 
computer implementation is small compared with the general 
form of sixth order digital notch-filter of thirteen 
coefficients. Therefore, one can implement the sixth 
order digital notch-filter (4-6) by direct canonic reali- 


zation as given by Fig. (4-42). 
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The corresponding state equations are given by: 


X4 (K) SEIT, (K il). + u (K) 

X1 (K*1) - x] (K) 

X» (K+1) = X1 (K*1) 

X (K+1) = xə (K*1) 

X 4 (K+1) = Xz (K*1) (4-7) 
Ze (Kr1) = X4 (K+1) 

Ze (K+1) = Xp (K+1) 

v (K) 5 X, (K+1) + x] (K) 


Oné Still TequUlrées the nOtEh frequency alt = OU 


0 
the equation (4-1) gives the sampling frequency T as follows 


T = se = سس شیم‎ = 1.388888889 x 10 ° sec. 
It is noted that the sampling time for sixth order is one 
half of that for the third order case. 

Figures (4-43) to (4-46) illustrate the time response 
of a sixth order digital notch-filter with input sinewave of 
60 Hz frequency. 

Figures (4-47) to (4-49) are frequency response plots 
of one thousand, two thousand and two thousand and five 
hundred iterations, respectively, corresponding coefficient 
a = 0.98298407. 

Itis noted that the sixth order digital notch- filter 
requires about fifteen thousand iterations to reach the 


steady-state. The steady-state gain frequency responses are 
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illustrated in Figs. (4-50) to (4-55) and some numerical 


values are tabulated in Table (4-4). 
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at 60 H2 Noch and 2500 iterations: 
a=0.98298407 
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One is lead to attempt to realize the sixth order 
30-7۲21 notch-filter by cascading three sections of second 


order as follows: 


u(K) V (KJ 
0131 Oly olf 
A ca E 


Figure (4-56) Realization of sixth order digital 


notch-filter by cascade. 


The steady-state gain frequency response is plotted 
in Fig. (4-57) and (4-58) for coefficient a=0.99 and a=0.995, 
respectively. 

Compare with the results of direct realization of 
Saben order notch-filter in Fics (452) and (4 53): the 
notch-gain in realization by cascading is deeper than in 
direct realization. But in passband it does not improve 
anything. It seems worse than direct form by introducing 
more noise; one can see gain frequency response behavior 
for coefficient a=0.9999999 between them in Figs. (4-59) 
and (4-60). 

It proves that the result is contrary with the 
statement of many current lectures: the pole-position of 
direct form is sensible than cascade form in realizing 


uen Order digital filter {Ref. 1]. The result in 





5105. (4-59) and (4-60) states that: the pole-position 
in cascade form is much more sensible than direct form in 
realizing sixth order digital notch-filter characterized 


by the transfer function (2-16). 
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Figure (4-57) Frequency response of sixth order 
notch-filter by cascade. 
a=0.99 
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Figure (4-59) Frequency response of sixth order 
nocchier by direct form. 
a=0.9999999 
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Figure (4-60) Frequency response of sixth order 


notch-filter by cascade form. 
a=0.9999999 
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4. Observation and Discussion 

In theory, by inspection, the magnitude response 
in equation (2-9) and the notch-width in equation (2-19) 
vn coefficient of m-th order digital notch-filter charac- 
terized by transfer function (2-16), one sees that when 
coefficient, a, approaches one, the magnitude of notch-gain 
and notch-width will be smaller [see Fig. (4-11)]. From 
the magnitude response plots, it is seen that theory and 
practice are in agreement whenever the steady-state output 
response is reached. This can be seen in the steady-state 
frequency responses in Figs. (4-12) to (4-17) of second 
order, Figs. (4-32) to (4-37) of third order and Figs. (4-50) 
to (4-55) of sixth order. Some of numerical results of 
steady-state response are tabulated in Tables (4-1), (4-2) 
and (4-3). 

For the case of intermediate state (e.g., state of 
one thousand iterations--the steady-state responses are not 
yet reached), the theory and practice appear to disagree. 
One can see this in frequency response plots in Figs. (4-2) 
to (4-10) of second order. These phenomena are a result 
of the random errors produced by truncating coefficients 
and round-off error in iteration computations in the simu- 
lating program due to finite register length of the computer. 
In this case, it is due to limitations of single precision 
arithmetic on IBM 360. A second explanation of the apparent 
deviation from the theoretical stems from error in computa- 


tion of notch-gain in computer program: 
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۱ _ Peak of input 
Gainag = 20 10210 (leak oF output) 


The peak of input is always one, but in the computer program, 
the peak of input and output are picked simultaneously and 

it is difficult to choose the sampling time such that it 
exactly corresponds to input peak; this is demonstrated in 


Fig. (4-56) 


time 





Figure (4 6 ٣ 


where the sampling time T is fixed for prescribed magnitude 
response of desired notch-filter. 
Suppose: 
06 1 ۶ 


the peak or output 222000002 


1909 





۱ 








he corresponding notch-gain will be: 
CAIN = - 75.06 0d 


Now let the peak close to one; suppose: 
the peak Of Input +0 
the peak of output - 0.0001 


GAIN * - 79.9] dB. 


The error thus realized approximately 6 dB for 
il) notceh-gain. -By the difficulty of picking exactly the 
peak of input as well as the peak of output, the actual 
notch-gains are certainly somewhat deeper than the values 
listed in Tables (4-1), (4-2), and (4-3). Also, the actual 
passband ripples should be smoother. 

It is known that the second order needs around 
six thousand iterations, the third order needs around ten 
thousand iterations and the sixth order needs around fifteen 
thousand iterations to the steady-state. It is seen by 
considering the time-respond plots of second, third and 
sixth order. It is obviously true since higher order in 
direct realization requires larger number of delays. 

It is noted that the passband-ripple is large 
around notch-frequency and smoother when going further. 

The results show that the higher order of digital 
notch-filter improves the passband-ripple (i.e., reduced 
ripple) and notch-gain [i.e., deeper notch. See Tables 


(4-2), (4-2), and (4-3)] and steady-state frequency response 
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& 


pisan Fics. (450) to (4-55) of second, third and 
sixth order, respectively. 

The sixth order notch-filter requires about four 
minutes to plot one steady-state gain frequency response 
using fifteen thousand iterations. 

One will see that the results in this chapter will 
give the easy way to work in optimum problems due to finite 


register length on computer in the next chapter. 


E CCD IMPLEMENTATIONS 

The charge-coupled device (CCD) was invented at the 
Bell Telephone Laboratory by Drs. Boyle and Smith in 1970. 
This new development has been receiving much attention 
by industry because of its important properties--low power 
consumption, low noise, and delay capability. Many success- 
ful applications have been investigated such as digital 
Signal processing especially in design of recursive digital 
filter. More details of CCD can be seen in many current 
lectures. The CCD hardware implementation of recursive 
digital notch-filter is beyond this research. The following 
is the sample of block diagram of second order CCD digital 


: 117 ۔ط 770 


*''Charge-coupled Device for Analog Signal Processing-- 
Recursive Filters Study," Engineer's Thesis by Varachai 
Imsa-ad, Naval Postgraduate School, December 1974. 
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Figure (4-61) Block diagram of a second order 


GGD digital  nəoteherfilter. 
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COEFETCİENT‏ گ۷ 


İn this chapter, a practical method is described for 
choosing the coefficients of a digital filter to meet 
arbitrary specifications of the magnitude characteristic. 

As shovn above, all special purpose computers have 
finite wordlengths which is a major source of noise in the 
implementation of digital filters. This is because of the 
necessity to truncation coefficient length and also because 
of rounding in the arithmetic operations, especially in 
the iterative multiplications. This kind of noise deter- 
mines the signal-to-noise ratio of the system. In some 
cases, the limitation of wordlength of coefficients may 
lead to a system which does not satisfy the original 
conditions. 

Especially, in the design of recursive digital notch- 
filters, one has to select the coefficients such that they 
satisfy the prescribed character o£ notch-gain, notchs 
vidth, and the smallest passband ripple. 

The m-th order digital notch-filter is characterized 
by the transfer function (2-16). It is seen that the 0 
gain, the notch-width, and the passband ripple are functions 
of a single parameter--the coefficient "a" between 0 and 1. 
It is possible to construct a performance function of 


single parameter in the optimization problem. 
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A. SEARCH PROCEDURE 

In the design. of an m th order digital notch filter; 
if one wishes the prescribed magnitude having the follow- 
ing characteristics: 

(i) Notch-gain is at least - 60 dB. 

(ii) Notch-width is at least one hertz. 
any coefficient, a, in )0,1( satisfies conditions (i) and 
(ii), and that value which produces the smoothest passband 
is the optimum value. 

In the method of searching, first, one selects the 
subinterval of )0,1( such that any coefficient belonging 
to this subinterval will satisfy the conditions (i) and 
(ii); second, searching is accomplished for the minimum 
error function of single variable for one of the 
subintervals. 

The subinterval can be found by computer experiment. 
Suppose for sixth order notch-filter one can consider 


Table (4-3); call (ај. аг) this subinterval: 


21 = 


OL pA SA hp hh hp اس کو ا‎ HSA AAA, 


As stated above, any coefficient, a, belonging to 
interval (as > a^) will satisfy the conditions (i) and (ii). 


For sixth order, the values of a, and a. are given as 


1 2 
[see Table (4-3)]: 


a, = 0.98 and a, = 0.995 


114 








For third order [see lable ld ol 


21 - 0.98 and a, - 0.99 


It is noted that the ideal (errorless) magnitude 
frequency response is equal to zero dB over the passband 
of the notch-filter. Due to the limited wordlength of 
computers, the actual magnitude frequency response will 


deviate from zero dB over the passband. If one defines: 
e (a, o.) = 1 - H (a, wj j (5-1) 


as the error function corresponding to coefficient, a, and 
radian frequency wj of sinusoidal input signal. The error 
mc tion € (a, ws) is defined in discrete domain which is 
considered as approximates the mean-squared error in con- 


tinuous domain of interval between 0 and We 


О). 
7 : | (a, نت‎ EO Je do (5-2) 
wj p co یں‎ 57: d 


€ (a; wj d 


wüere H. (a, Ws) is ideal response and has magnitude equal 
to one. 


Rewrite (5-2) as 


Y. 
1 


es | 2 : 
€ (a, ws) H x 1 H (a, w, ) | do; (5 3) 


the system is stable, e (a, ws) has to be convergent. 


1 





Let o. Scan over the frequency range of the passband of 
netch-filter, the error performance or the total error 


function, is defined as 


M nu 
E(a) = = ce” (a, o.) 
ızı 
ör 
S 2 
Ea (aera (5-4) 
әд 
where aj < a < a, and N is aa finite integer number. 


The optimization problem turns out to be the least 
square method of single parameter. The optimum coefficient, 
a, 1s one that makes the error performance E(a) minimum. 

The following algorithms are described to search for the 
minimum of Efa). First, the interval (ај. a) is subdivided 


into N subintervals such that 


K(a,-a,) 
EM o 


)0- ی۰ی‎ ии 
d +یا‎ a 1 7: 


K 1 N 
where K = 1, 2, ..., N. 
Then compare the values of E (ax) and E (441) Corresponding 


to successive value of ay and ak4]? where 


K(az-a.) (5:3) 
2 ٦ 25 
acı “alı N vu am m” 
pr B(ay) 5 E (24,1) 
the next step will compare E (ak+1) and E (Ax 47) and store 


coefficient a+]: 
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8 


0" 07۶0 (ay) — (3,4) , the next step will compare 


E (ay) and E (a and store coefficient a the execution 


K*2) K’ 
continues until the last subinterval of (ay > a^) is reached. 
The final result will be printed out as the optimum 
coefficient corresponding to the smallest value of error 

pemrormance E (a) for aj > 8 € a». 

The greater the value of N, the higher accuracy one 
can get but the price is paid by spending more computation 
time. 

The plots of error performance are illustrated in Fig. 
(5-1) of second, third, and sixth order digital notch- 
filters. One Sees that the error performance is monotoni- 
cally decreasing with increasing coefficient, i.e., the 
greater the coefficient, the smoother the response in 
passband. This result can be extended for many order of 
digital notch-filter characterized by the transfer function 
(2-16). The optimum coefficient now is so easy to pick and 
it is the coefficient a,. 


2 
for Sixth order a 


0.995 


for third order 2 0.99. 


p DESCUSSTƏN AND CONCLUSION 

One sees that the error performance of Single parameter 
is unimodal. The optimum problem of this kind of function 
has been accomplished by many methods well-known such as 
Fibonacci Search, Golden Section, etc. The Fibonacci Search 


is known as the most efficient in optimization technique of 














sixth Order 
E 5. Third order 


us ae Second order 


750 Cli ۴ 
¿00 UNITS INCH. 


Plot Of total deviation of second! 

third and sixth order:notch 11067 

between a=0.5 and a=0.999, 

The region in circle is exaggerated 
Jn Fig. (522), 


1: 
y 





Figure (5-1) 
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unimodal objective function of single parameter. The method 
used here--referred to by some authors as the uniform or 
factorial search--compares formatly with the Fibonacci 
Search in that it takes more computation time than, but in 
this case of error-performance E(a), in the author's 
opinion, is the most efficient one. In Spite of the cost 

of computation time, the method is simple and the result 

1s precise. 

The results carried out in this chapter agree with 
computer experiments in Chapter IV by considering the 
numerical value in Tables (4-1), (4-2), and (4-3). Higher 
orders of notch-filter caused deeper notch-gains and smoother 
response in passband. It is noted that the smoother pass- 
band is correct for coefficient close to one. It is seen 
in Fig. (5-2) for noise-performance in the region close 
tO One. | 

The value of the optimum coefficient, a, increases with 
the order of filter provided the conditions (i) and (ii) 
are fixed. The higher order of notch-filter, the closer a 
approaches one, the constraint is placed on a as mentioned 
in section D of Chapter II, due to the finite wordlength, 
e.g., the single precision on IBM 360/70 has 32-bit 


wordlength. The inequality (2-23) gives 


l - a, = À > 2791 


Or 


a, « 1 - 277İ . 0.999999999 5 
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Figure (5-2) Noise performance in the region 
of coefficient close to one: 
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The value (1 - 2°°!) is the upper limit of coefficient, a, 
in theory. In practical view, the upper limit is smaller 
than (1 - تک‎ by the reason of error introduced in 
arithmetic operation and error in representation as stated 
above. Besides, in the implementation of notch-filter, 
the optimum coefficient corresponding the larger wordlength 
in register of computer will be closer to unity than 
the shorter one. For example, in implementation of sixth 
order of single precision on IBM 360, the upper limit of 
coefficient is 0.999, for coefficient exceeds this limit, 
the system begins to degenerate. It is the same limit for 
second and third order. One can see this phenomenon 
numerically in Tables (4-1), (4-2), and (4-3) and the 
behavior of changing in steady-state frequency response 
lots of second, third and sixthiorder in previous Chapter: 
One sees that for a = 0.9999, the notch-gain less 
than 3 dB in magnitude and the passband ripple is unacceptable. 
It is foolish to accept them as the notch- filters: 
Furthermore, ۶۲0۰۰43 > — the system may become 
the: oscillator due to error in rounding since it could 
make the poles of their transfer functions exactly on the 


unit Circle Gl -z plane:- 
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VI. CONCLUSION 


The filter design was used to design the digital notch- 
alters, -The effect of order ol Tilter ompnocch san) 
notch-width, and passband ripple was investigated. 

From results of this study, the simplest configuration 
of digital notch-filters was found; hence implementation 
will be economical in hardware. The notch-frequency, the 
notch-gain and the notch-width are very easily controllable 
by varying only one coefficient in realization as well as 
the agreement between the predicted and practice sampling 
time. This property is usable in an array of notch-filters 
to eliminate interference signals of variety of frequencies. 

It was seen that with selected coefficient the pass- 
band is very smooth (for sixth order less than 0.05 dB) 
and the notch-width is very narrow (less than one hertz). 
The improvement can be achieved by increasing the order of 
filters provided the coefficient close enough to the optimum 
coefficient which was found by Search procedure. The opti- 
mum computer program was used and the result was precise. 
This program required large computer time. That is the 


reason some authors do not like it. 
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